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Abstract 

Some intriging connections between the properties of nonlinear 
■ Q noise driven systems and the nonlinear dynamics of a particular set 

of Hamilton's equation are discussed. A large class of Fokker-Planck 
O Equations, like the Schrodingcr equation, can exhibit a transition in 

their spectral statistics as a coupling parameter is varied. This tran- 
T^j sition is connected to the transition to non-integrability in the Hamil- 

ton's equations. 
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In this paper we will be concerned with diffusion processes on W 1 de- 
© scribed by the set of coupled stochastic differential equations 

O dq t (t) = -d l $(c L )dt+^dW' t (t), i = l,...,n, (1) 

where 3>(q) is a potential bounded from below, the W l (t) are uncorrelated 
Wiener processes, and g is a diffusion coefficient. In this case the evolution 

© of the probability density p(q, t) on 5ft n , is described by the Fokker-Planck 

l_j equation 

g %>=§Ap + V>V<&). (2) 

Using the time separation ansatz p(q, t) = p(q)e~ xt / g , we can write Eq. (2) 
as an eigenvalue equation Cp\(q) = —\p\(q), where C = ^-A + gV $ + 
<?V<I> • V. After the change of basis p(q) = e _ */ 9 ^(q) we obtain 

W* A (q) = A* A (q), (3) 

where Ti = —e®^ 9 C e~®^ 9 = — ^-A + $(q), is a Hermitian Schrodinger type 
operator with the transformed potential l> = ^(V^) 2 — § V 2 ^. The problem 
of solving Eq. (2) has been reduced to the problem defined by equation (3). 
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For small g the WKB solutions of equation (3) are given by 



* A (q) = £c a |VS Q r^exp i-S a {d,\)\ , (4) 

where the S a (q, A) are the solutions of the Hamilton-Jacobi equation ^(VS a ) 2 + 

= A. The solutions of this equation are given by the integrals S* a (q, A) = 
/ q p a • dq where the integration in along the classical trajectories of Hamil- 
ton's equations of motion 

dH . dH 
P = -W q= 5p' (5) 

with 

#(p,q) = ^P 2 + ^(q)- (6) 

The time reversal symmetry of equations (5) with Hamiltonian (6) insures 
that the eigenfunction (4) are real since the solutions of the Hamilton-Jacobi 
equation will come in pairs ±S a . The dynamics of (5) determine the solution 
of (3) through (4), and the solutions of (2) are given by 

p A (q,t)=exp(- $ ^)* A (q). (7) 

Thus the properties of the Fokker-Planck equation (2) are connected to 
the dynamics of the system with Hamiltonian (6) in a manner somewhat 
analogous to the relation of a quantum mechanical system to its classical 
counterpart. 

One question we might ask is how the behavior of (2) is affected by the 
degree of chaos in the equations of motion (6). Such effects, in the quan- 
tum mechanical case ((5) affecting (3)), are often referred to as quantum 
chaos, which is usually defined as the characteristics of quantum systems 
whose classical analogues exhibit chaos. The statistical properties of the 
eigenvalues of such systems are such characteristics, and the level spacing 
distribution P(S), giving the probability of level separation S (measured 
in units of the local mean spacing), provides one such statistical property. 
Berry & Tabor [1] have shown that nearly all quantum systems whose clas- 
sical analogues are integrable will have a Poisson level spacing distribution 
P(S) = exp(— S), indicating the statistical independence of neighboring en- 
ergy levels. On the other hand, it is now understood that the eigenvalues of 
systems whose classical analogues are chaotic exhibit level repulsion. That 
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is, P(S) — > as S — > 0.[2] It is expected that systems with time- reversal sym- 
metry whose classical analogues are globally chaotic will have a Wigner level 
spacing distribution, P(S) = ttS/2 exp(— 7r5 2 /4), indicating a linear level re- 
pulsion as S — ► 0.[3] Since the eigenvalues of the Fokker-Planck operator, C, 
with potential $ are the negative of the eigenvalues of a Hamiltonian, 7i, 
with potential $, the spectral statistics of the Fokker-Planck equation (2) 
would then be expected to provide a signature of the dynamics of the equa- 
tions of motion (5). To explore these ideas Millonas and Reichl[4] studied a 
family of two-dimensional Fokker-Planck equations with potentials 

*e(x,y) = 2x A + ^y 4 + exy(x-y) 2 . (8) 

The system needed to be at least two-dimensional in order to observe chaos 
in equations (5). When e = the system is completely integrable, since it 
decouples into two one-dimensional systems. They observed the transition 
(as e is varied) in the level spacing statistics of the Fokker-Planck operator 
as the dynamics of equations (5) changes from completely integrable (e = 0) 
to almost globally chaotic (e = 0.14). Stochastic chaos can then be defined, 
at least for the case of diffusion in a time-independent potential, as the 
properties of stochastic systems described by Eq. (2) when the equations of 
motion ( 5) exhibit chaos. In particular, given a family of potentials $ e where 
the dynamics of (5) varies from globally integrable to globally chaotic as e 
is increased, we would expect the spectral spacing distribution of the A's 
to exhibit a corresponding transition from Poisson to Wigner level spacing 
statistics. 

An entirely separate problem is the question of the direct physical rel- 
evance of the dynamics of (5) to the underlying microscopic dynamics as 
described by (1). One thing is clear: chaos in (5) is emphatically not related 
to chaos in the dynamics generated by (1) with g = 0, that is q = — V<l?(q). 
When there is no noise the individual trajectories just follow the gradient 
of the potential along the route of steepest descent stopping at any local 
minimum in $, so what would normally be considered the underlying mi- 
croscopic dynamics is trivial, and never chaotic. Thus, there is no simple 
physical relationship between the dynamics of (5) and the dynamics of (1). 
A deeper analysis shows that eqs. (5) are the imaginary-time equations of 
motion for the most probable, or optimal trajectories. These ideas can be 
extended to the case where there is no detailed balance, but in that case no 
meaningful analytic continuation of the most probable trajectories is pos- 
sible. There are still however the optimal trajectories which obey a set of 
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Hamilton's equations with Hamiltonian Ti = l/2(p — A) 2 + 1/2V • A, where 
A is the nongradient force field which replaces — <9*<£(q) in Eq. 1 in the 
case where there is no detailed balance. The long-time behavior of systems 
with or without detailed balance can be calculated in the low-noise asymp- 
totic limit from a knowledge of the optimal trajectories with energy A = 0. 
These are the instanton trajectories which line on the unstable manifold 
of the Hamiltonian system. This manifold is smooth as a consequence of 
the center manifold theorem, so chaos will not play a role. However this 
manifold may have singular projections onto the configuration space in the 
case where detailed balance is broken, resulting in a rich nonlinear behavior 
of nonequilibrium stationary states. [5] The most surprising result presented 
here is than even in the case where there is detailed balance, and the under- 
lying dynamics is completly integrable, chaos will play a role in determining 
the time-dependent properties of such systems. It appears that there is a 
deep analogy between quantum dynamics and stochastic dynamics through 
their relationship to the properties of these conservative dynamical systems. 
This connection, once made, opens up the study of stochastic processes to 
a whole range of new tools and concepts. 
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